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ABSTRACT 

Wind-blown bubbles, from those around massive O and Wolf-Rayet stars, to super- 
bubbles around OB associations and galactic winds in starburst galaxies, have a dom- 
inant role in determining the structure of the Interstellar Medium. X-ray observations 
of these bubbles are particularly important as most of their volume is taken up with 
hot gas, 10 5 < T(K) < 10 8 . 

However, it is difficult to compare these X-ray observations, usually analysed in 
terms of single or two temperature spectral model fits, with theoretical models, as 
' real bubbles do not have such simple temperature distributions. Spectral fits, and the 

\Q , properties inferred from them, will depend in a complex way on the true temperature 

■ distribution and the characteristics and limitations of the X-ray observatory used. 

C*~) ' In this introduction to a series of papers detailing the observable X-ray properties 

of wind-blown bubbles, we describe our method with which we aim to solve this 
problem, analysing a simulation of a wind-blown bubble around a massive star. 

Our model is of a wind of constant mass and energy injection rate, blowing into a 
uniform ISM, from which we calculate X-ray spectra as would be seen by the ROSAT 
Qh PSPC. Analysing these spectra in the same way as a real observation would be, we 

compare the properties of the bubble as would be inferred from the ROSAT data with 
^ ■ the true properties of the bubble in the simulation. 

^ | We find standard spectral models yield inferred properties that deviate signifi- 

ed , cantly from the true properties, even though the spectral fits are statistically accept- 

able, and give no indication that they do not represent to true spectral distribution. 
For example, single temperature spectral fits give best fit metal abundances only 4% 
>^ . of the true value. A cool bubble has best fit temperatures significantly higher than 

$H ' a bubble twice as hot. These results suggest that in any case where the true source 

spectrum does not come from a simple single or two temperature distribution the "ob- 
served" properties cannot naively be used to infer the true properties. In this situation, 
to compare X-ray observations with theory it is necessary to calculate the observable 
X-ray properties of the model. 

Key words: Methods: data analysis - Methods: numerical - ISM: bubbles - X-rays: 
interstellar 
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1 INTRODUCTION 

Bubbles blown in the Interstellar Medium (ISM) by mas- 
sive stars are a common astrophysical phenomenon. X-ray 
observations can provide information regarding the density, 
metal abundance, temperature, ionisation state and phys- 
ical structure in the hot bubbles surrounding Wolf-Rayet 
(WR) and O stars (Wrigge, Wendker & Wisotski 1994), Lu- 
minous Blue Variables (LBV's) such as r\ Carinae (Corco- 
ran et al. 1995) and planetary nebulae (PN) (Kreysing et 
al. 1992; Leahy, Zhang & Kwok 1994; Arnaud, Borkowski 



& Harrington 1996; Leahy et al. 1996). On the larger scale, 
superbubbles are created by the sum of the winds and SN 
within OB associations (Belloni & Mereghetti 1994) and gi- 
ant star forming regions in young starburst galaxies. Su- 
perbubbles within starburst galaxies such as M82 eventu- 
ally break out the galaxy to form spectacular galactic winds 
(Watson, Stanger & Griffiths 1984; Heckman, Armus & Mi- 
ley 1987). In many cases the X-ray emission probes different 
regions of the object in question to that revealed by optical 
observations, increasing the importance of the X-ray data. 

Analytic solutions to the development and structure of 
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astrophysical bubbles must rely on simplifying assumptions, 
and increasingly attention has turned to the use of numeri- 
cal hydrodynamics. These simulations have been enlighten- 
ing with respect to the nonlinear processes occurring, with 
some degree of quantitative agreement with observation, but 
generally lack predictive power. This is partially due to the 
difficulty in comparing them with observations, in particular 
X-ray observations. 

The problem is that X-ray observations are usually 
analysed by fitting a single or two temperature spectral 
model to the observed spectra (see for example the refer- 
ences above), and the best-fit results are used to infer the 
physical properties of the object. 

However, for wind-blown bubbles such as those men- 
tioned above, the true situation is more complex, and the 
results of the spectral fits may be influenced by, for exam- 
ple: projection of different physical regions along the line of 
sight; the presence of a wide range of temperatures; interstel- 
lar absorption; unknown or non-standard elemental abun- 
dances; non-ionisation equilibrium conditions; low numbers 
of observed photons and the limitations of the current X-ray 
telescope optics and detectors. 

All of these make the interpretation of what is normally 
a one or two temperature spectral fit to the data difficult to 
relate to the underlying physical conditions, and conversely, 
make it difficult to predict the observable properties of a 
model or simulation. 

To our knowledge there has been no study of the in- 
fluence of the complexities mentioned above on the best-fit 
properties of a spectral fit to the observable X-ray data, 
and in particular not for wind blown bubbles. As we shall 
show, the combination of the physical effects above and the 
properties (and limitations) of real X-ray observatories, can 
significantly affect the results of simple spectral fitting. 

Previous authors (Weaver et al. f 977; Zhekov & 
Perinotto f996) have calculated theoretical X-ray spectra 
from their f-D models, but did not consider particular in- 
struments or fit models to those spectra. In general only 
X-ray luminosities are calculated {e.g. Volk & Kwok 1985; 
Mellema & Frank 1995; Garcia-Segura & Mac Low 1995). 

The aim of this paper is to introduce a method of 
analysing numerical simulations in the same way as actual 
X-ray observations are analysed, i.e. predict the observable 
X-ray properties. This method can be applied to a wide 
range of phenomenon where X-rays are important, from 
wind-blown bubbles around WR stars and PN, through the 
larger bubbles around clusters of massive stars to starburst- 
driven galactic winds. 

We simulate a wind blown bubble using a 2-dimensional 
hydrodynamic code, concentrating on the properties of the 
hot X-ray emitting gas. The hydrodynamic model is used 
to generate artificial X-ray spectra and images, in particular 
simulated ROSAT spectra. We then analyse these spectra in 
the same way as real ROSAT spectra would be, in order to 
determine what the observationally determined properties 
of the bubble would be, and how those relate to its actual 
structure. This synthesis is necessary to a) determine the 
physical processes that are observationally important, and 
b) allow a direct comparison between observation and theory. 

Our model of a wind-blown bubble is deliberately cho- 
sen to be the simplest applicable model with an analytic 
solution, in order to simplify the analysis of our results, and 



avoiding added complications that a more physically accu- 
rate model of a wind blown bubble (e.g. Garcia-Segura, Mac 
Low & Langer 1996) would introduce into the interpreta- 
tion of our results. Later papers will consider more realistic 
models, with additional physics such as time varying energy 
and mass injection rates, along with spatial variation of the 
X-ray properties. This will be necessary to understand the 
properties of the extended emission from galactic winds (see 
for example Strickland, Stevens and Ponman 1997). 

In Section |^ we describe the numerical code used to 
produce the results shown in Section []. Section ^ discusses 
the implications of these results, and we briefly sum up in 
Section 5l 



2 NUMERICAL METHOD 

The 2-D numerical simulations presented in this paper 
were performed using Virginia Hydrodynamics- 1 (VH- 
1), written by John Blondin and co-workers (Blondin et 
al. 1990). VH-1 is based on the piecewise parabolic method 
(PPM) of Colella & Woodward (1984), which is a third order 
accurate finite difference algorithm. For a more detailed de- 
scription of the code see Stevens, Blondin & Pollock (1992). 

Radiative cooling is handled assuming optically thin gas 
with solar abundance ratios and ionisation equilibrium. The 
cooling curve for the temperature range 4.0 < logT < 8.6 
was generated using the Raymond-Smith plasma code (Ray- 
mond & Smith 1977), assuming solar element abundances. 
Because cooling times at T ~ 10 K are short, we place 
a limit of the numerical timestep to prevent more than a 
2% energy loss per computational cell per timestep. For 
T < 10 4 K the cooling is set to zero. 

For the simulations described below VH-1 was run in 
cylindrical coordinates (r, z), assuming symmetry around 
the z-axis. In practice this gives rise to numerical artifacts 
along the coordinate axes, but these are small and do not 
affect our results. The computational grid consisted of 400 x 
400 cells of equal size, spanning a physical distance of r max = 
Zmax = 2.0 x 10 19 cm. 

2.1 Simple wind-blown bubble model 

As a simple realistic model of a wind-blown bubble we con- 
sider the analytical model of Castor, McCray & Weaver 
(1975) and Weaver et al. (1977). A wind of constant termi- 
nal velocity % and mechanical luminosity tw = |Mw"w 
(where Mw is the mass loss rate of the star) blows into 
a medium of uniform density po- This generates a bubble 
with a structure divided into three distinct regions: an outer 
shock, separating the undisturbed ambient ISM from a shell 
of swept up and shocked ISM; a bubble of hot shocked stel- 
lar wind material, bounded at its outer edge by a contact 
discontinuity between it an the shell of shocked ISM, and at 
its inner edge by a reverse shock between it and the third 
innermost region of the freely expanding stellar wind. After 
several thousand years the shell of swept up ISM cools and 
collapses down to form denser, cold (T ~ 10 4 K) shell, which 
is the source of the the optical emission associated with the 
nebula. 

In the standard Weaver et al. (1977) model thermal con- 
duction leads to evaporation of matter off the dense shell of 
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Figure 1. Logarithm of the gas number density during the simulation at t = 3500, 7700, 10170 and 14630 yr. At t = 3500 yr the bubble 
has suffered no significant radiative energy loss. Shell collapse is underway at 7700 yr, approximately the time of maximum soft X-ray 
luminosity, and has just finished at 10170 yr. The bubble then enters the self-similar phase, its properties at t = 14630 yr being typical 
of this stage. 



swept-up ISM into the hot bubble interior, cooling it and 
significantly increasing its density. In common with many 
of the hydrodynamical simulations of bubbles we do not in- 
clude conduction. Magnetic fields may suppress conduction, 
even at very low B-field strengths that otherwise are dy- 
namically unimportant (Soker 1994; Band & Liang 1988). 
Given the limited observational knowledge on the state of 
magnetic fields and conduction within bubbles we choose to 
ignore both! 

We shall consider two simulations, both with L\y = 
6.3 x 10 37 ergs _1 . To obtain different temperatures in the 
hot bubble we vary the wind mass loss rate between sim- 
ulations. The low mass loss rate simulation has Mw = 
5 x 10 -5 Mq yr -1 . In the high mass loss rate simulation 



Mw = 1 x 10 -4 Mq yr -1 , giving a bubble of half the temper- 
ature of the low mass loss rate simulation. Mass and energy 
are added to cells within r = 3 x 10 18 cm at each timestep. 
The ambient medium is assumed to be uniform and totally 
ionised, with a total number density of no = 10 cm -3 . 



2.2 Simulated X-ray data 

We shall only consider simulated ROSAT PSPC (Posi- 
tion Sensitive Proportional Counter) data. Although the 
gas proportional counter's spectral resolution of AE/E m 
0.43(£/0.93)-°- 5 (FWHM, with E measured in keV) is low 
compared to a mission such as ASCA, wind-blown bubbles 
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Figure 2. As Fig. [l]but greyscalcs of log 10 T(K). 



Table 1. Model parameters for the two bubble simulations 



Parameter 
Wind luminosity 
Mass loss rate M\y 
Wind velocity u\y 
ISM ambient density no 
Distance to bubble 
Absorbing column TVh 
Grid size 
Cell size 



Value 

6.3 X lO^crgs" 1 

5 X 1CT 5 or 1 X 10- 4 M Q yr- 1 

2000 or 1414 km s" 1 

10 cm -3 

2kpc 

3.16 X 10 21 cm" 2 

2. 10 19 cm X 2. 10 19 cm 

5. 10 16 cm X 5. 10 16 cm 



are soft X-ray sources and R OS AT has more sensitivity than 
ASCA at low energies. 

For a detailed discussion of the ROSAT satellite see 
'The ROSAT user's handbook' (Briel et al. 1994). 



3 RESULTS 

As the low and high mass loss rate simulations only differ in 
the density and temperature of the shocked wind, we shall 
concentrate on describing the low mass loss rate simulation 
below. Section B.4 describes how the results of the high mass 
loss rate bubble differ from those given below. 



3.1 Bubble growth 

There are three definable stages of bubble growth seen in 
the simulation. These are i) before the shell cools, ii) during 
shell cooling and collapse, and iii) self-similar growth after 
shell collapse with a thin cold (T ~ 10 4 K) shell. The 1- 
dimensional analytic solutions for the first and last of these 
stages are presented in detail in Castor et al. (1975) and 
Weaver et al. (1977). 
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Figure 4. Radial profiles of density, temperature, velocity and pressure (solid lines) from the low mass loss rate simulation, compared 
to the predictions from of the Weaver et al. (1977) model (with conduction: dashed line, no conduction: dotted line), at t = 14630 yr. 
Note that the Weaver model only applies for R > Ri, the radius of the wind termination shock. 



Initially the swept up ISM is shock-heated to 5.5 < 
log T(K) < 6.0. The shell is thick (see Figs. |l| and |), and a 
strong emitter of extreme Ultraviolet (EUV) radiation and 
soft X-rays, as can be seen from the 0.1-2.4 keV luminos- 
ity (Fig. |§]). The major coolant is radiation in the UV-EUV 
rather than X-rays, the UV-EUV luminosity being of order 
a magnitude greater than the soft X-ray luminosity before 
shell collapse. The luminosity rises rapidly with time, as the 
bubble sweeps up and heats more ISM. The rate of increase 
of luminosity decreases after t ~ 4000 yr, the X-ray lumi- 
nosity peaking at Lx = 1.9 x 10 36 ergs _1 at t w 6300 yr 
as the shell begins to cool. The UV-EUV luminosity peaks 
later, at « 8100 yr as the shell cools further out of the X-ray 
band and becomes denser. The peak UV-EUV luminosity of 
9.1 x 10 37 ergs -1 briefly exceeds the wind energy input. Af- 



ter shell collapse X-ray luminosities are approximately two 
orders of magnitude below the UV-EUV luminosity, both 
remaining essentially constant for the duration of the simu- 
lation. 

Following shell collapse shell densities are typically sev- 
eral hundred to a few thousand particles per cubic centime- 
tre, with T « 10 4 K. 

In the absence of heat conduction and evaporation off 
the cool shell, we would simplistically expect the bubble 
interior to have a uniform low density and high temperature. 
In this case the shocked- wind material has the same pressure 
as predicted by Castor et al. (1975) and Weaver et al. (1977), 
but the temperature determined by the reverse shock (the 
termination shock of the freely expanding wind) . 

In practice, the density rises in as the shell is ap- 
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Figure 3. Bubble X-ray luminosity as a function of time in the 
ROS/IT 0.1-2.4keV band. 




10 



0.01 0.1 
T (keV) 



10 



Figure 5. The emission measure as a function of temperature 
{EM T = J nj? x dV T ) for the entire bubble at t = 15000 yr. 



gations seen on the inside surface of the shell, shear motions 
arise between the faster bubble interior and the shell, leading 
to swirling motions along the interface and the introduction 
of cooler, denser material into the hot bubble. 

In a perfectly spherically symmetric bubble, the lack 
of a tangential velocity component between the faster- 
expanding bubble interior and the shell would prevent such 
stripping of material off the shell. In our simulations, the 
bubble-shell interface is corrugated by instabilities from 
early on in the simulation, presenting faces not totally per- 
pendicular to the flow in the bubble interior, and leading to 
mixing. 



3.1.1 Shell instabilities 

The instabilities of the shell seen in Figs. |and| have impor- 
tant consequences as they lead to the introduction of cooler, 
denser material into the bubble interior, hence modifying 
the X-ray emitting properties of the bubble. 

The shell should be stable against Ray leigh- Taylor in- 
stabilities, as it is constantly decelerating, suggesting that 
the instabilities are Vishniac instabilities (Vishniac 1983). 

The initial seed perturbation is numerical artifact, aris- 
ing when the forward shock first appears at the start of the 
simulation, and is due to the orthogonal nature of the com- 
putational grid and the finite size of the energy injection 
region. 

Does this instability, and the cooler material it injects 
into the hot shocked wind region, render the results of these 
simulations untenable? In real system such as wind-blown 
bubbles around individual stars or superbubbles, both the 
wind energy injection rate and the ISM will be far from uni- 
form, so instabilities and inhomogeneities in the shell are 
to be expected. Note that conduction would introduce far 
more cool material into the bubble interior than these in- 
stabilities do (Fig. Q). Mass loading of the flow by ablation 
of clumps and cloudlets (c.f. Hartquist et al. 1986) will also 
add cool dense material. We are therefor not worried that 
our 2D simulations do not exactly match the ID analytic 
prediction, as the mass introduced into the bubble interior 
is by no means unphysically large. 



3.2 Bubble structure and properties at 
t = 15000 yr 

For a more detailed study of the observable X-ray proper- 
ties we shall investigate the bubble at t = 14630 yr (here- 
after 15000 yr for convenience). The bubble at this time has 
settled down into self-similar expansion with a thin shell, 
the period described by the standard solution of Castor et 
al. (1975) and Weaver et al. (1977). 



proached, and the temperature drops, although not to the 
extent expected for true conduction. This can be seen in 
both the 2D images of Figs. ^-|^ and the radial profiles of 
Fig. ^. This is due to eddies and swirling motions along the 
shell-bubble interface mixing material from the dense shell 
into the hot bubble interior. The outward velocity in the 
shocked- wind is higher than the velocity at which the shell 
expands into the ISM (Fig. [|), and coupled with the corru- 



3.2.1 Comparison with the Weaver et al. model 

The properties of our simulated bubble at t — 15000 yr and 
the predictions of the Weaver et al. (1977) for the same in- 
put parameters are encouragingly similar, if we modify the 
standard Weaver et al. model to account for the lack of ther- 
mal conduction in the VH-1 code (Fig. ^ and Tables ^ and 

To illustrate the effect that conduction would have on 
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Table 2. Bubble properties at t = 15000 yr. The X-ray fluxes 
in the ROSAT 0.1-2.4keV band and ROSAT PSPC count rates 
are given for zero column (intrinsic) and Nn ■■ 
(absorbed). 



Parameter 


Value 




Reverse shock radius R± 


5.3 X 10 18 cm 


Forward shock radius R2 


1.17 X 10 19 cm 


Total energy radiated 


1.057 x 10 49 


erg 


Thermal energy 


1.169 x 10 49 


erg 


Kinetic energy 


6.824 x 10 48 


erg 


Shell velocity 


126-166 km s 


-1 


Total volume 


6.702 x 10 57 


cm 3 


Total mass 


34.4 M© 




Total emission measure 


4.409 x 10 60 


cm -3 


Luminosity (0.1-2.4 keV) 


6.55 x 10 34 ergs -1 


fx (int.) 


1.36 x 10- 10 


erg s 


PSPC count rate 


35.3 photons 


s" 1 


fx (abs.) 


3.54 x 10~ 12 


ergs - 


PSPC count rate 


0.3156 photons s 



Table 3. Bubble properties at t = 
Weaver et al. (1977) model. 

Parameter 

Reverse shock radius Ri 
Forward shock radius R2 
Thermal energy in shocked wind 
Kinetic energy in shell 
Shell velocity 

Total volume (including shell) 
Total mass (including shell) 
Luminosity" (0.1-2.4 keV) 
Luminosity 6 (0.1-2.4 keV) 

a Including thermal conduction. 
b Without conduction. 



15000 yr as predicted by the 



Value 

4.2 X 10 18 ( 



5.67 x 10 48 erg 
143 km s- 1 
5.6 X 10 57 cm 3 
28.7 M e 

2.39 X 10 35 ergs -1 
8.94 X 10 32 ergs -1 



the structure of the bubble, the properties of the standard 
Weaver et al. bubble with conduction (with thermal con- 
ductivity of the form K{T) = CT 5/2 , where G = 1.2 X 
10 -6 erg s -1 cm" 

-i K -7/2 ( Spitzer 1962 )) 

are also given in 

Fig. I and Table |. 

Note that the analytical predictions only apply for R > 
Ri, the radius of the reverse shock, although they are shown 
in Fig. ^extending all the way to the center of the bubble. In 
Castor et al. (1975) and Weaver (1977) the free-wind region 
inside Ri is ignored, which as we shall discuss below, does 
lead to minor differences between simulation and theory. 

Due to the way we have injected mass and energy onto 
the computational grid, the region inside Ri corresponds to 
the solution of Chevalier and Clegg (1985), rather than a 
point-like source of mass and energy. 

The Weaver et al. model's shell is much thinner than the 
azimuthally averaged profile from the simulation, as the shell 
is corrugated and hence spread out over radius. Inspecting 
the 2-D data the peak densities in the simulations agree well 
with the density predicted by the analytical model, imply- 
ing that the dense shell is not significantly unresolved. The 
bubble in the simulation is slightly larger than predicted, 
but this is not surprising as this reflects the "head-start" 
we have given the simulated bubble by injecting mass and 
energy over a finite input region. 

Densities (temperatures) in the shocked-wind rise 
(drop) near the shell, compared to the analytic prediction, 



due to the mixing of coo ler material into the bubble as dis- 
cussed in Section. |3.1.l[ This increases the bubble X-ray 
luminosity from L x = 8.94 x 10 32 ergs -1 in Table. || for a 
non-conductive bubble to an appreciable 6.55 x 10 34 ergs -1 , 
only a factor ~ 3 less than for a bubble with thermal con- 
duction. 

The thermal pressure in the shocked wind in our simu- 
lation is less than that predicted by the Weaver et al. model. 
In the analytic solution it is assumed that the energy con- 
tained in the free-wind region is negligible, and that all the 
energy in the shocked wind is thermal. In practice this is 
not the case, although the majority of the energy is thermal 
energy in the shocked wind. The velocities in the shocked 
wind are appreciable (see Fig.^). Coupled with the appre- 
ciable energy loss that accompanied shell collapse this leads 
to lower than expected thermal energy in the shocked wind, 
and hence lower pressure. 

3.2.2 Simulated R OS A T spectra 

We generate and fit theoretical X-ray spectra by the fol- 
lowing process. The 2-Dimensional density and temperature 
data from VH-1 are rotated around the axis of symmetry 
to produce a temporary 3-D dataset. The MEKAL (Mewe, 
Kaastra & Liedahl 1995) plasma code is used to calculate 
the X-ray emission from each volume of gas, assuming ioni- 
sation equilibrium and solar element abundances. It is also 
assumed the gas is optically thin to its own radiation, and 
that any absorption the X-rays suffer is extrinsic to the 
source. External absorption is modeled as a uniform fore- 
ground screen, using the coefficients given in Morrison & 
McCammon (1983). 

The raw high resolution X-ray spectrum generated by 
the plasma code (see Fig. |^) is corrected for a typical ab- 
sorbing column of 3.16 x 10 21 cm -2 , appropriate for the 
wind-blown bubble NGC 6888 (Wrigge, Wendker & Wisot- 
ski 1994), and then folded through the spectral response of 
the ROSAT PSPC, assuming a source distance of 2kpc for 
the bubble. 

No X-ray background is included, effectively assuming 
a perfect background subtraction. We estimate that the soft 
X-ray background would have a PSPC count rate of « 0.19 
counts s _1 over the area of sky occupied by the bubble at 
t = 15000 yr, based on the ROSAT PSPC observation of 
M82 (Strickland et al. 1997). The detector count rate due to 
the bubble itself is 0.32 counts s -1 for an absorbing column 
of 3.16 x 10 21 cm -2 . 

Note that the time to reach ionisation equilibrium, 
ticq ~ 10 12 n -1 s (c.f. Masai 1994) is greater than the bub- 
ble's age over much of the volume. As our aim is to see 
whether fits to ROSAT spectra can accurately reflect the 
true bubble properties, the assumption of ionisation equlib- 
rium is not a problem. 



3.2.3 A simple estimator of the fitted temperature. 

Are there easier ways of predicting the observable properties 
of a model than generating and fitting artificial X-ray spec- 
tra? This method is also instrument specific, as it depends 
on the spectral characteristics of the instrument used, and 
so predicted results for ROSAT are not directly transferable 
to, e.g. ASCA. 
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If we assume we can correctly estimate the absorbing 
column and metallicity, either from a spectral fit or from 
observations at other wavelengths, then estimating a char- 
acteristic temperature is all we need do. 

We might expect the fitted temperature to be some flux- 
weighted average of the true temperature range, and the 
emission measure roughly the total emission over the range 
of temperatures ROSAT is sensitive to. 

Over the energy range 0.005-15.0 keV the intrinsic flux 
(i.e. unabsorbed) weighted average temperature < Tew >= 
(4.9^'g) x 10~ 3 keV, i.e. UV emission from very cool mate- 
rial at the shell-bubble interior interface. In practice all the 
mid and extreme UV emission is absorbed in the intervening 
ISM. In addition the instrument response further limits the 
"observed" temperature. 

If we assume that we can accurately remove the ef- 
fects of interstellar absorption on the X-rays (i.e. fitting 
the column correctly) then to construct a predicted ROSAT 
temperature we need only weight each temperature by the 
count rate in the ROSAT PSPC due to each tempera- 
ture component, ignoring absorption. This temperature is 
< Trqsat >= 0.09881°;^ keV. Given the true tempera- 



ture distribution in Fig. g the ROSAT emission-weighted 
average temperature is more a reflection of the instrument 
response than a fair estimation of the temperature of the 
emitted radiation or the general state of the bubble. 

< Trosat > is instrument specific, defeating our object 
of trying to produce an instrument-independent estimator 
for the temperature. This may well be impossible, in such 
a situation as this where the temperature distribution can 
not be described by a single, characteristic temperature. 

The emission measure should follow from the normalisa- 
tion of the spectrum, and will depend on the best-fit values 
for the other parameters which determine the shape of the 
model spectrum. For the purposes of comparison between 
the results of the spectral fitting below and the "true" val- 
ues, the total emission measure between 0.01 < T( keV) < 
10.0 is 3.77 x 10 58 cm -3 . Between 0.1 < T(keV) < 2.4, it is 



3.2.4 Spectral fitting 

For an assumed exposure time of 3000 s we obtain a sim- 
ulated ROSAT PSPC spectrum containing ~ 1000 counts. 
To assess the effect of real photon statistics, and the con- 
sequent variation in best-fit parameters, on our simulated 
spectra we generated ten Poisson realizations, fitting each 
using the Starlink X-ray analysis package Asterix. The 
number of counts per bin is sufficient to allow the use of \ 2 
fitting. 

Three different spectral models were considered: a stan- 
dard single temperature absorbed hot plasma model; a two 
temperature model, with both components having the same 
absorbing column and metal abundance, and finally a dif- 
ferential emission measure model. 

Single temperature models are widely used to charac- 
terise X-ray emission, and it is only when there are sufficient 
counts to show that a single temperature model is a bad fit 
that more complex models are used. It is therefor sensible to 
fit such a model, despite knowing that the true temperature 
distribution of the X-ray emitting plasma is far from being 
single temperature (see Fig. ^). 



A two temperature model would be the next level of 
complexity, naively a soft component for cooler denser gas 
near the shell and a hot component for the bubble interior. 

The differential emission model attempts to incorporate 
the emission of gas at a wider range of temperatures by spec- 
ifying the emission measure EM = J ri^dV is a power law 
in temperature, i.e. EM oc T 7 , between two temperatures 
Tiow and Thigh. 

We quote results for models with all parameters fitted 
for (including the column and the metallicity) and for fits 
with the metallicity fixed at solar abundance. Although in 
practice we know the absorbing column applied, fits with 
it fixed were statistically unacceptable. For the purpose of 
display and interpretation, the best-fit results for each of the 
ten Poisson realisations were averaged. The quoted errors 
are the statistical deviations of the best-fit results, and not 
the fitted confidence regions. Typically, the 68% confidence 
regions (for one parameter of interest) were smaller or of 
order the deviations quoted. 

The results for the fits to the integrated bubble spec- 
trum are given in Tables. ^, |E] and ^j. 

The best fit models look very similar when convolved 
with PSPC instrument response (i.e. as they are when com- 
pared with the data), despite the very different fit parame- 
ters (Fig. Q). The only noticeable differences are the width 
of the peak at channels equivalent to E ~ 0.8 keV, and how 
the model fits to points between 1.5 < E(keV) < 2.0. 

The reason the single temperature, metallicity fixed at 
solar, spectral fits are bad is that they fail to fit these higher 
energy points at all. Freeing the metallicity to fit allows 
the shape of the spectrum to be changed, the tempera- 
ture increases to better fit the higher energy points, and the 
metallicity drops to change the relative shape of the spec- 
trum around the Iron-L complex (effectively the only lines 
ROSAT spectral fitting is sensitive too). Similarly the ab- 
sorption column changes the shape of the spectrum. In the 
fits with metallicity free to fit, the increased temperature 
move the spectrum towards higher energies, removing the 
requirement to absorb as much low energy flux, hence low 
best fit columns. 

The two temperature and differential emission measure 
models give better fits, as the hot component (or hot compo- 
nents in the case of the differential emission measure model) 
fit the high energy points present in the spectrum, free- 
ing the cool component to fit at a lower temperature. The 
available trade-off between the various components in these 
multi-component models make the best-fit parameter less 
constrained, and insensitive to the metal abundance. 

We shall defer further interpretation and discussion of 
these results, and what they would lead us to infer about 
the bubble's properties, to Section. 4.1.5. 



3.3 Spectral variation over time 

The spectral properties of the bubble after shell collapse will 
be different to those prior to shell cooling. At early times, as 
discussed above, the shell is the major source of X-rays, and 
dominates the emission. It is only after the shell cools from 
its initial temperature of T ~ 10 6 K to the final tempera- 
ture of T ~ 10 4 K that the X-ray emission from the bubble 
interior, i.e. the shocked wind itself, becomes significant. 
This is graphically illustrated by Figs. ^]andj|. Note how 
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o X— ray spectrum at t = 14630 yr 

O p 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — i 




0.5 1 1.5 2 

E (keV) 

Figure 6. X-ray spectra for the bubble at t = 15000 yr, un- 
convolved with any instrument response. The solid and dashed 
lines correspond to absorbing columns of A/h = 0.0 and 3.16 X 
10 21 cm -2 respectively. For comparison an unabsorbed T = 
0.1 keV spectrum of arbitrary normalisation is shown dotted. Note 
that despite being soft the bubble's spectrum does have a hard 
tail. 



Table 4. Average result of fits to simulated ROSAT spectra using 
a single temperature MEKAL hot plasma model with absorption. 
The errors quoted are the deviations between the best fit results. 
See text for details. Note the significant improvement in fit when 
metal abundance is free to fit. 



Parameter 


Single temperature fits 


Units 










N H 


7.8 ±0.4 


10 21 cm" 


2 








EM 


0.018 ±0.005 


10 60 enr 


3 








T 


0.163 ±0.006 


keV 










Z 


1.0 (Fixed) 


Zq 










fx (obs.) 


2.86 ± 0.07 


10~ 12 erj 


;s" 


i 


cm 


-2 


fx (int.) 


7.18 ± 1.84 


10" 10 erj 


;s _ 


i 


cm 


-2 


X 2 (d.o.f) 


39 ±5.6 ( 19 ) 












N H 


1.9 ±0.4 


10 21 cm" 


2 








EM 


0.0019 ±0.0012 


10 60 enr 


3 








T 


0.45 ± 0.07 


keV 










Z 


0.04 ± 0.02 


Zq 










fx (obs.) 


3.29 ± 0.10 


10" 12 er{ 


;s" 


i 


cm 


-2 


fx (int.) 


0.17 ±0.07 


10" 10 erj 


;s _ 


i 


cm 


-2 


X 2 (d.o.f) 


18.5 ± 4.8 ( 18 ) 













after shell collapse the averaged temperature contributing to 
the photons detected by ROSAT, Trosat, increases again. 
This is the emission from shocked wind itself, previously of 
too low a level compared to the X-ray emission from the 
shell to be noticeable. 

Fig. [] shows the best fit temperatures from single tem- 
perature fits to a set of simulated ROSAT spectra over time, 
for both metal abundance fixed at solar and abundances 



Channel space ROSAT PSPC spectrum 
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Figure 7. Channel space ROSAT PSPC spectrum of the bubble 
at t = 15000 yr shown with the three average best fit models 
(single temperature, abundance fixed at solar - solid line; single 
temperature, abundances fitted for — dashed; two temperature, 
abundances fixed - dotted) The two temperature and differential 
emission measure models give very similar ROSAT spectra. The 
spectrum has no noise, the error bars represent the size of typical 
Poisson uncertainties for an exposure time of 3000 s. 



Table 5. As Table ^ but fitting with a two temperature MEKAL 
model. 



Parameter 


Two temperature fits 


Units 






JVh 


3.5 ±2.8 


10 21 cm" 2 






EMi 


0.04 ±0.10 


10 60 cm" 3 






Ti 


0.23 ±0.08 


keV 






EM 2 


0.00013 ± 0.00006 


10 60 cm" 3 






T 2 


1.28 ±0.30 


keV 






Z 


1.0 (Fixed) 


Zq 






fx (obs.) 


3.32 ±0.10 


10 12 ergs 1 


cm 


-2 


fx (int.) 


13.7 ±32.9 


10~ 10 ergs" 1 


cm" 


-2 


X 2 (d.o.f) 


12.0 ±4.6 ( 17 ) 








N H 


3.7 ±2.9 


10 21 cm" 2 






EM! 


0.013 ±0.031 


10 60 cm" 3 






Ti 


0.24 ±0.11 


keV 






EM 2 


0.00005 ± 0.00006 


10 60 cm~ 3 






T 2 


1.45 ±0.37 


keV 






Z 


6.7 ±5.0 


Zq 






fx (obs.) 


3.31 ±0.11 


10 -12 ergs -1 


cm" 


-2 


fx (int.) 


33.3 ±80.9 


10-10 ergs" 1 


cm" 


-2 


X 2 (d.o.f) 


11.5 ±4.1 ( 16 ) 









free to fit. As generating 10 Poisson realisations for each of 
~ 50 spectra, fitting and finally averaging the results would 
be prohibitively time consuming, we have fitted the spectra 
with no Poisson noise applied. Comparison of this with the 
more rigorous method used earlier shows that this gener- 
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Figure 8. The flux weighted average temperature Tew of the 
bubble (see Sect. |l.l| ), as a function of the age of the bubble, in the 
band 0.005-15.0 keV. At early times the emission is predominantly 
from the swept-up and shock heated ISM. After this shell cools 
the total bubble emission is still dominated by the shell, but the 
radiation is now UV rather than soft X-rays. Dashed lines show 
the r.m.s. deviation of Tew 

Trosat compared to best— fit temperatures 
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Figure 9. The ROSAT flux-weighted aver age t emperature 
Trosat as a function of bubble age (see Sect. 3.2.3), compared 
the temperatures obtained from single temperature fits to the 
PSPC spectra (Diamonds: fits with metal abundance fixed at so- 
lar, crosses: fits with metal abundance free to fit). Error bars are 
68% confidence in one parameter of interest. The dashed lines 
show the r.m.s. deviation of Trqsat- 



Table 6. As Table ^ but fitting with a differential emission 
measure model, where the emission measure at any temperature 
i?M( T j oc T 7 , between two cut-off temperatures Ti ow and Thigh- 
This model attempts to recreate the multi-phase origin of the 
emission. 



Parameter 


Diff. emission measure fits 


Units 








JV H 


3.6 ± 2.3 


10 21 cm -2 








EM 


0.00376 ± 0.00524 


10 60 cm -3 








7 


-0.23 ±0.36 










^low 


0.0118 ± 0.0345 


keV 








Thigh 


7.4 ± 2.2 


keV 








Z 


1.0 (Fixed) 


Z Q 








fx (obs.) 


3.37 ± 0.13 


in— 12 - 

10 iz ergs 


l 


cm" 


_2 


fx (int.) 


1.46 ± 2.03 


10-10 ergs" 


1 


cm" 


-2 


X 2 (d.o.f) 


15.2 ±6.3 ( 17 ) 










JVh 


3.0 ±1.8 


10 21 cm" 2 








EM 


0.00264 ± 0.00287 


10 60 cm" 3 








7 


-0.29 ±0.39 










Tow 


0.0071 ± 0.0193 


keV 








Thigh 


8.7 ±4.6 


keV 








z 


1.16 ± 1.53 


Zq 








fx (obs.) 


3.38 ±0.12 


10 12 ergs 


1 


cm" 


-2 


fx (int.) 


0.81 ± 1.33 


10 — 10 erg s~ 


1 


cm" 


-2 


X 2 (d.o.f) 


14.2 ± 5.7 ( 16 ) 











ates essentially identical results, although the values of \ 2 
obtained are misleadingly low. 

As can be seen from Fig. |^ the best fit temperatures 
fail to reflect the temperature changes we know are occur- 
ring. Intriguingly, the best fit metallicity does show a sys- 
tematic trend between t ~ 5000 and 10000 yr. Two temper- 
ature fits, able to fit the shape of the spectrum better, do 
show an initial drop in both temperatures to a minimum at 
t ~ 5000 yr, followed by a increase to a constant value after 
t ~ 9000 yr (although the hotter component only levels out 
after t ~ 15000 yr). 



3.4 The high mass loss rate simulation 

To further investigate how spectral fitting to the observed 
X-ray data depends on the the properties of the bubble, we 



3.2 



have repeated the detailed analysis described in Section 
on the higher wind mass loss rate simulation. 

Increasing Mw to lO~ 4 Af0yr _1 results in a bubble 
with the same size at t = 15000 yr, identical cold shell prop- 
erties, but a shocked wind with a density (temperature) a 
factor 2 higher (lower). 

Single temperature models give best-fit results very sim- 
ilar to those given in Table ^ for the low mass loss rate sim- 
ulation, except the best fit temperatures are significantly 
higher in this case: T = 0.41 ± 0.09 keV (metal abundance 
fixed) or T = 0.60 ± 0.13 keV (metal abundance fitted for). 
This is despite the average temperature within the shocked 
wind being half that of the low mass loss rate simulation. 
Again absorption columns (and when fitted for, metallici- 
ties) deviate systematically from the true values. 

The two temperature spectral models give best fits very 
similar to the results given in Table H the best fit tempera- 
tures being Ti = 0.26 ± 0.07 keV and T 2 = 1.47 ± 0.31 keV. 
The only real differences are a) the relative fraction of the 
total emission measure in the hot component has increased 
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(i.e. the spectrum appears to be harder), and b) there is 
evidence for two fit minima, one with low Nh and emission 
measure, the other a high Nh with high emission measure. 



4 DISCUSSION 

4.1 Properties inferred from the spectral fitting 

We shall concentrate on the results from the low mass loss 
rate simulation, mentioning only where the spectral fit re- 
sults from the high mass loss rate simulation differ. 



4-1.1 Single temperature models 

Single temperature spectral fits give misleading information 
on the state of the bubble, not surprising given the multi- 
component temperature structure and the detector charac- 
teristics, nor are they sensitive to the true spectral changes 
that occur during the bubble's growth (see Fig. This 
suggests that the best-fit values are strongly affected by the 
multicomponent structure and the instrument response. 

The confidence levels from the fits misrepresent the true 
uncertainties in the fitted parameters. In particular the fit- 
ted temperatures are apparently well constrained, although 
the emission comes from a very broad temperature distri- 
bution. The spectral fits give little clue to the true tem- 
perature distribution, as the fits with metallicity free to fit 
have acceptable % 2 ■ A naive hope that in a situation such 
as this, with a broad, non-peaked temperature distribution, 
single temperature fits would show poorly constrained best- 
fit temperatures is not justified by these results. 

Absorbing columns and metal abundances are both sys- 
tematically mis-fitted. It is clear that by fitting a multicom- 
ponent spectrum with a single temperature model we intro- 
duce a major systematic effect. Single temperature fits with 
metal abundance fixed at solar (Table ^|) are generally unac- 
ceptable, with reduced \ 2 ~ 2, perhaps indicating that the 
model is not a good representation of the data. The hydro- 
gen columns deviate significantly from the true absorbing 
column. 

Freeing the metal abundance to fit does give good fits 
to the data, but at the expense of giving best-fit parameters 
that bear little relation to the bubble's true properties. The 
best fit metal abundances are less than one twentieth so- 
lar, and apparently strongly constrained. This is interesting 
given the current debate over the accuracy of X-ray deter- 
mined abundances (c.f. Bauer & Bregman 1996). 

The fl ux-we ighted average temperature Trosat derived 
in Section 3.2.3 is not a good estimator of the best-fit tem- 
perature, as can be seen in Fig. ^. An average temperature, 
even sensibly weighted, does not reflect how spectral fitting 
works, and will not give acceptable results. 



4-1.2 Two temperature models 

These (Table ||) give significantly better fits to the data 
than the single temperature models and are slightly better 
than the differential emission measure models. The absorb- 
ing columns fit closer to the true value. The fits are insen- 
sitive to the metallicity, and would not force us to believe 



they are significantly different from solar, unlike the single 
temperature models. 

That the best-fit temperatures in the two temperature 
model fits to the low and the high mass loss rate simula- 
tions are so similar, with only the relative normalisation of 
the two components varying, suggests the temperatures are 
determined more by the ROSAT PSPC's response than the 
true temperature distribution of the source. 



4-1.3 Differential emission measure models 

From Fig. ^| it is clear that for T < 0.7 keV the emission 
measure is well represented by a power law in T, with a slope 
of 7 » -1.8. Between T = 0.1 keV and T = 0.7keV the true 
emission measure falls by nearly two orders of magnitude, so 
a differential emission measure model with this slope should 
provide a good approximation to the spectrum. 

The results of the model fits (Table ^|) are therefor sur- 
prising, given the best fit values for 7 ~ —0.25! The com- 
bination of material hotter than T ~ 1 keV and the energy 
dependent response of the ROSAT PSPC must bias the fit 
significantly. The slope is clearly wrong within the energy 
range 0. 1-2.4 keV ROSAT is sensitive to, let alone extended 
to the cut-off temperatures claimed by the fit. 

4-1-4 The effect of longer exposure times 

Are any of the effects above due to low numbers of photons? 
Repeating the above analysis with assumed exposure times 
of 10000 s, giving simulated spectra with > 3000 counts, 
gives results very similar to those in Section [| The confi- 
dence regions quoted on the individual fits are smaller, but 
the average best fit result is consistent with those quoted 
above, and are statistically acceptable fits. 

It therefor appears that the systematic deviations of the 
spectral fit results from the true values is due to fitting an 
intrinsically complex spectrum with a simplistic model, and 
not due to poor photon statistics. 

4-1.5 Inferred properties 

Can we infer any of the true bubble properties from the 
spectral fits? In particular, the true luminosity, densities, 
pressures and thermal energy are interesting quantities that 
we would like to know in addition to the temperature if 
X-ray observations are to be of any use in understanding 
the object. 

The 0.1-2.4 keV fluxes predicted from all the models are 
reasonably accurate before correction for absorption (com- 
pare the real values from Table |^ with the results from the 
spectral fitting, Tables^ - ^). The single temperature models 
get the intrinsic luminosity very wrong, due to the incorrect 
absorption column, either overestimating it by a factor ~ 5 
(metallicity fixed at solar) or underestimating by an order 
of magnitude in the case of the fits with the metallicity free 
to fit. On average the two temperature models overestimate 
the intrinsic flux by an order of magnitude, although with 
a large uncertainty due to the large variation in best-fit col- 
umn. Similarly the differential emission measure models also 
have a large scatter in inferred intrinsic luminosity. 

Assuming we know the volume of the emitting plasma 
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Figure 10. An X-ray image of the bubble at t = 15000 yr, in the 
0.1-2.4keV energy band. The image has been smoothed with a 
Gaussian of FWHM 25" to approximate the ROSAT PSPC PSF. 
Units are counts s _1 arcmin - 2 . 



V, the root mean square electron density is n c = 
yf\EM/V), where EM is the (volume) emission measure 
emission measure obtained from the spectral fit. Then the 
thermal pressure and energy are P w 2n e kT and Eth ~ 
3n c kTV. 

Given the X-ray surface brightness (Fig. |l(J) we can 
estimate the volume. The inferred properties are relatively 
insensitive to V as P and n c are proportional to V _1//2 and 
Eth oc V 1/2 . 

The bubble is clearly limb-brightened (Fig. |lo| ), the 
emission coming predominantly from a thin shell near the 
shocked- wind/cold shell interface. From a radial profile of 
the ROSAT surface brightness, and correcting for the point 
spread function, the inferred radial thickness of this shell is 
~ 20", the PSF-corrected volume Korr = 9.8 x 10 56 cm 3 . 
By way of comparison, the total volume within the bubble 
inferred from the X-ray image is Vbub = 6.0 x 10 57 cm 3 . Note 
that this is less than the volume quoted in Table ^ as we are 
not including the cold shell. 

The density, pressure and thermal energy inferred from 
the single temperature models, using Vcorr and Vbub, are 
given in Table m 

For the two temperature models, it would be sensible 
to associate the cool component (with a high emission mea- 
sure) of the spectral fit, to cool, dense material near the 
shocked- wind/cold shell interface, given the observed limb 
brightening. The hot component in the fit, with a total emis- 
sion measure two to three orders of magnitude less than the 
cool component, could well represent emission from the hot 
rarefied gas in the bubble interior. Assigning volumes Vcorr 
and Vbub to the cool and hot components respectively we ob- 
tain densities, pressures and thermal energy content given 
in Table fjj 

From the simulation we know the true properties of all 
the material in the bubble. The density in the plasma con- 
tributing the majority of the 0.1-2.4 keV flux (all the ma- 



terial with volume emissivities within an order of magni- 
tude of the maximum emissivity) varies between 2.5 < n c < 
28.2 cm -3 , averaging 7.8 cm -3 . In the hotter bubble interior 
the electron density is n c » 0.09 cm" 3 . The shocked wind 
and cold shell are practically isobaric at P = 1.3 x 10 -9 dyn 
cm -2 , and the majority of the total thermal energy in the 
bubble of -Eth = 1.17 x 10 49 erg is in the hot shocked wind 
material that occupies most of the bubble volume. 

Comparing the inferred properties of the bubble from 
Table |7| with the true properties (see for example Fig. ^) we 
find: 

• The single temperature spectral fits, assuming the emis- 
sion come from a shell of volume V CO rr, give, on average, 
estimates of the density and pressure (in the plasma that 
dominates the X-ray emission detected) that are within a 
factor 2 of the true values. 

• Thermal energies inferred from single temperature 
models generally underestimate the true thermal energy in 
the bubble. The single temperature models are dominated 
by the cool gas, occupying only a fraction of the total vol- 
ume, whose properties are not typical of the bubble as a 
whole, and miss the hotter gas that contains most of the 
thermal energy. 

• Two temperature models give a better idea of the bub- 
ble properties, although the temperature of the hot compo- 
nent is significantly lower than the true temperature within 
the shocked wind. 

• By attempting to fit a intrinsically multicomponent 
spectrum with simplistic single temperature model, the best 
fit hydrogen column and metal abundances can be signifi- 
cantly in error. In more complex models this effect is re- 
duced, but not eradicated. For the three spectral models 
considered the fitted temperatures give little indication of 
the true temperature distribution. 

• Spectral fits give a good estimate of the observed energy 
flux, but extrapolating the intrinsic source flux is error prone 
due to its strong dependence on the assumed absorption 
column, and may be an order of magnitude out. 

• The results of the spectral fitting give little indication 
of the true variation in properties within the bubble, i.e. 
the multiple temperature components. In particular, com- 
ponents with temperatures outside the primary bandpass of 
the instrument are not represented by the spectral fitting. 
The differential emission measure model fails to give a bet- 
ter fit than a two temperature model, although on paper it 
should perform better. There would be little reason to be- 
lieve that the two temperature model was in fact an over 
simplification of the true situation. 

• The systematic error introduced by under-modeling the 
true spectrum will dominate over any statistical error on 
the best-fit parameters. The best fit may be a statistically 
good fit, but does not give any indication that it is a poor 
representation of the true properties. 

It is interesting to note that despite the difference be- 
tween the true bubble properties and those for a bubble with 
thermal conduction (compare the radial profiles in Fig. ^), 
the properties inferred from the ROSAT PSPC two temper- 
ature spectral fitting are generally consistent with a bubble 
with conduction: 

• The predicted luminosity of a conductive bubble is 
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Table 7. Bubble properties as would be inferred from the ROSAT 
spectral fits to the low mass loss rate simulation. The volumes 
refer to the volumes defined in the text. 



Model 


Volume 


n c 

( cm -3 ) 


P 

(dyn 


cm 


- 2 ) 


(erg) 




IT, Z = 1 


Vcorr 


4.3 


2.2 X 


10" 


-9 


3.3 x 


10 48 




Vbub 


1.7 


9.0 x 


10" 


■10 


8.2 x 


10 48 


IT, Z free 


Vcorr 


1.4 


2.0 x 


io- 


9 


2.9 x 


10 48 




Vbub 


0.6 


8.0 x 


io- 


■10 


7.3 x 


10 48 


2T, Z = 1 
















(cool) 


Vcorr 


6.6 


4.8 x 


10" 


9 


7.1 x 


10 48 


(hot) 


Vbub 


0.15 


6.1 x 


io- 


■10 


5.5 x 


1Q 48 


2T, Z free 
















(cool) 


Vcorr 


3.1 


2.9 x 


io- 


9 


4.3 x 


10 48 


(hot) 


H>ub 


0.09 


4.3 x 


10" 


■10 


3.9 x 


10 48 



4-1.7 Summary of "observed" properties 

Fitted absorption columns can vary wildly in a case with a 
complex spectrum, biasing estimated emission measures are 
luminosities. It might be safer not to fit for the column, and 
use estimates from optical or Hi measurements. 

The metallicity in single temperature models is badly 
underestimated. Given the current debate over X-ray deter- 
mined abundances (Bauer & Bregman 1996), where much 
of the attention has concentrated on the hot plasma codes 
rather than the effect of multi-temperature gas, this is a 
very interesting result. Is this a peculiarity of this particular 
model, or a general property? Does this also affect higher 
spectral resolution instruments such as ASCA or the future 
XMM7 

As single temperature fits fail to represent to true emit- 
ted spectrum they generally underestimate the thermal en- 
ergy in the bubble, as this is primarily in hot, tenuous gas 
with much lower emissivity than the cool gas. We note that 
this could provide an explanation for the apparent discrep- 
ancy reported by Magnier et al. (1996) between the thermal 
energy in the superbubble N44 as derived from ROSAT and 
ASCA spectral fits, where the X-ray derived thermal en- 
ergy is a factor 6-10 lower than expected on the basis of the 
Weaver et al. model. However, we have not demonstrated 
that our results apply over such a wide range of parame- 
ter space, and the alternative explanation that the stellar 
wind energy input may have been overestimated may well 
be true. We are continuing this work to specifically model 
superbubbles in future papers. 

Nevertheless, it is safe to say that the results of simple 
spectral fits to what are in reality more complex sources de- 
pend strongly on both the spectral response of the detector 
and the spectral distribution of the source. 



slightly higher than that of the simulated bubble, but the 
spectral fits generally overestimate the intrinsic luminosity 
of the bubble. 

• The temperatures predicted by the model with thermal 
conduction are T ~ 1 keV in the shocked wind, dropping to 
T ~ 0.2 keV in the denser material near the boundary with 
the cold shell (Fig. ^J), very similar to the best fit tempera- 
tures. 

• The mass flux off the cold shell driven by thermal con- 
duction naturally leads to a region of warm X-ray emitting 
gas of the right density just inside the cold shell. 

• Only the density of the hot component is inconsistent 
with the value predicted by the model with thermal conduc- 
tion. The inferred density in the hotter material is too low 
by a factor ~ 5. 



4-1-6 High mass loss rate simulation 

The true properties of the bubble in the high mass loss rate 
simulation are practically identical to those of the low mass 
loss rate bubble. The only differences are that the density 
and temperature in the shocked wind are a factor of two 
different. The intrinsic X-ray luminosity in the 0. 1-2.4 keV 
band is 3% higher, but most importantly the ROSAT count 
rate is double that of the low mass rate simulation. 

The intrinsic luminosity is again dominated by the cool, 
dense material near that shocked wind/cold shell interface, 
whose properties are almost identical between the two sim- 
ulations. 

The count rate depends on both the intrinsic spectrum 
and the detector's spectral response. In this case the cooler 
bubble interior brings the emission from the shocked wind 
more within the ROSAT spectral response, giving more high 
energy photons in the PSPC spectra, and hence making 
them seem harder. This explains the higher best-fit tem- 
peratures of this cooler bubble, compared to the low mass 
loss rate simulation. 

This is further evidence, if needed, that it is necessary to 
predict the observable properties of a model when comparing 
observation with theoretical models. 



4.2 Limitations and assumptions 

We have not included several physical effects that are po- 
tentially important in real stellar wind blown bubbles: ther- 
mal conduction, magnetic fields, time-dependent ionisation, 
non-solar abundance ratios, mass loading and interaction 
with the winds from the main sequence and red giant stages 
of the star. Mass and energy were injected onto the com- 
putational grid in a manner more suitable for superbubbles 
around OB associations and in starburst galaxies. For the 
spectral fitting, we have ignored the soft X-ray background. 

These omissions do not invalidate or reduce the sig- 
nificance of the results contained in this paper. Our aim 
is to illustrate the need for modeling the X-ray emission, 
and analysing in as similar a manner to real observations 
as possible. In this introductory paper we have avoided 
the potential complications that additional physics would 
introduce into the interpretation of the results, and take 
the known properties from relatively simple model, cal- 
culate what would be seen in a real observatory such as 
ROSAT, and ask "can we infer the true properties from the 
X-ray emission resulting from those properties?" Relaxing 
the above assumptions only make the need for this method 
of direct comparison more urgent. 
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4.3 Future work 

We intend to extend this work, relaxing the assumptions 
of constant mass and energy injection, and a homogeneous 
ISM, to model the properties of superbubbles and galactic 
winds in starburst galaxies. As a general method and phi- 
losophy it can, of course, be applied to a wider range of 
astrophysical bubbles. 

In future we shall consider the spatial variation of the 
X-ray properties, as done for example for the ROSAT ob- 
servation of M82's galactic wind by Strickland, Ponman & 
Stevens (1997). 

We shall also extend this work to other X-ray obser- 
vatories such as ASCA or XMM, to make use of the supe- 
rior spectral resolution of these instruments. For the higher 
signal-to-noise spectra that XMM will make available, this 
method of direct comparison will be absolutely necessary! 



5 CONCLUSIONS 

We have shown that in order to compare X-ray observations 
to theory, it is necessary to consider the observable X-ray 
properties of the theory. The results of a spectral fit are a 
complex function of the the density and temperature dis- 
tributions of the source, absorption, the properties of the 
detector used and the spectral fitting procedure. As such 
they should not be considered as "real" values, but as char- 
acteristic values, and specific to the instrument used. The 
normal method of fitting a simplistic model to the observed 
data, and then treating the best-fit parameters as the real 
properties can easily give answers an order of magnitude out 
from the truth. 

This technique will allow the first direct comparison be- 
tween observation and theoretical models of superbubbles 
and starburst driven outflows. 

We would like to thank Trevor Ponman, Robin Williams 
and the referee for constructive criticism. DKS and IRS ac- 
knowledge financial support from PPARC. This work was 
performed on the Birmingham node of the Starlink net- 
work. 
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